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Localized precipitation and runoff on Mars 

Edwin S. Kite, 1,2 Timothy I. Michaels, 3 Scot Rafkin, 3 Michael Manga, 1,2 and 
William E. Dietrich. 1 

Abstract. We use the Mars Regional Atmospheric Modeling System (MRAMS) to sim- 
ulate lake storms on Mars, rinding that intense localized precipitation will occur for lake 
size >10 3 km 2 . Mars has a low-density atmosphere, so deep convection can be triggered 
by small amounts of latent heat release. In our reference simulation, the buoyant plume 
lifts vapor above condensation level, forming a 20km-high optically-thick cloud. Ice grains 
grow to 200 /im radius and fall near (or in) the lake at mean rates up to 1.5 mm/hr wa- 
ter equivalent (maximum rates up to 6 mm/hr water equivalent). Because atmospheric 
temperatures outside the surface layer are always well below 273K, supersaturation and 
condensation begin at low altitudes above lakes on Mars. In contrast to Earth lake-effect 
storms, lake storms on Mars involve continuous precipitation, and their vertical veloc- 
ities and plume heights exceed those of tropical thunderstorms on Earth. For lake sizes 
10 2 - 5 - 10 3,5 km, plume vertical velocity scales linearly with lake area. Convection does 
not reach above the planetary boundary layer for lakes <C10 3 km 2 or for atmospheric 
pressure > O(10 2 ) mbar. Instead, vapor is advected downwind with little cloud forma- 
tion. Precipitation occurs as snow, and the daytime radiative forcing at the land sur- 
face due to plume vapor and storm clouds is too small to melt snow directly (< +10 
W/m 2 ). However, if orbital conditions are favorable, then the snow may be seasonally 
unstable to melting and produce runoff to form channels. We calculate the probability 
of melting by running thermal models over all possible orbital conditions and weight- 
ing their outcomes by probabilities given by long-term integrations of the chaotic dif- 
fusion of solar system orbital elements [Laskar et al., 2004]. With this approach, we de- 
termine that for an equatorial vapor source, sunlight 15% fainter than at present, and 
snowpack with albedo 0.28 (0.35), melting may occur with 4% (0.1%) probability. This 
rises to 56% (12%) if the ancient greenhouse effect was modestly (6K) greater than to- 
day. 



1. Introduction 

Evidence for runoff on Mars shows it to be patchy in both 
space and time [Kraal et al, 2008; Williams, 2007; Weitz et 
al., 2008; Fassett & Head, 2008a; Hynek et al, 2010; Carr & 
Malin, 2000], so perhaps past precipitation was also patchy. 
Because patchy surface vapor sources cannot persist in equi- 
librium with a dry atmosphere [Richardson & Soto, 2008; 
Soto & Richardson, 2008], vapor would have to be supplied 
from an environment not in equilibrium with surface con- 
ditions. Such environments can be transient, such as an 
impact lake, or long-lived, such as the base of a wet-based 
ice-sheet. They can be high-temperature, such as fumaroles 
(or a lava flow advancing over snowpack), or involve only 
moderate temperatures, such as groundwater discharge. 

Here we use a mesoscale model to explore the atmospheric 
response to one example of a non-equilibrium vapor source: 
an ephemeral lake on a cold desert Mars. We track the fate 
of vapor supplied by the lake from release, through cloud 
formation, to precipitation, and consider whether the re- 
sulting snow will melt and provide runoff to form channels. 
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Lake size, solar luminosity, lake geometry, and atmospheric 
pressure all affect the results. Only idealized results are pre- 
sented: a companion paper {Kite et al. [submitted]; hence- 
forth Paper 2) uses the same model for a case study of the 
Juventae plateau inverted channel networks [Weitz et al., 
2008]. 

Low volumetric heat capacity makes the Mars atmo- 
sphere's response to lake vapor release similar to tropical 
moist convection on Earth, so we borrow ideas from trop- 
ical meteorology to understand our results (e.g. Emanuel 
[1994]). Figure 1 shows the low-pressure lake effect: con- 
densation of a small amount of vapor in a thin atmosphere 
can produce strong convection, which in Earth's thick at- 
mosphere would require condensation of a large amount of 
vapor and correspondingly high water surface temperatures. 

Localized precipitation on a cold desert planet is nor- 
mally transient precipitation. A warm, wet patch connected 
to the global atmosphere will lose water to cold traps else- 
where on the planet [Richardson & Soto, 2008]; the water 
table will withdraw to the subsurface because of evapora- 
tive losses [Soto & Richardson, 2008]. In the absence of an 
external heat source, evaporative and radiative cooling will 
cause any lake to quickly freeze [Lorenz et al., 2005; Conway 
et al., 2010]. (These arguments do not apply to springs, nor 
proglacial discharge of subglacial meltwater. In these cases, 
under cold conditions, any given parcel of water will freeze 
over, but a sustained vapor source can nevertheless exist at 
the discharge site.) As the ice thickens, ice surface temper- 
ature will fall and the lower saturation vapor pressure will 
cause the rate of vapor release to greatly decrease. For real- 
istic external heat sources, the lake lifetime is still short. For 
example, consider an impact-generated lake near the freez- 
ing point overlying shocked basalt that is initially at 1000 
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°C. The lake is assumed to be well-mixed by waves driven 
by lake-effect and impact-thermal storms, and convection 
driven by bottom heating. Icing-over is inevitable when the 
heat flow from the interior of the lake toward the surface is 
less than the evaporative and radiative losses at the surface. 
The maximum time before icing over, t, is therefore 

_ D(T b - T)c b p b 

~ (EL vap + oT*) 1 ' 

where D is the depth of pervasive fracturing within the 
rock cjecta, T b = 1273K is the initial temperature of the 
basalt, T = 278K is lake surface temperature, c b = 840 
J/kg/K the specific heat capacity of the basalt, p b — 2000 
kg/m 3 the density of the fractured basalt, E = 2 kg/m 2 /hr 
the evaporation rate, L vap = 2.5 x 10 J/kg the latent heat of 
vaporisation, and a = 5.67 x 10 -8 W/m 2 /K 4 is the Stefan- 
Boltzmann constant. For D — 100m we obtain ~4 Earth 
years: a geological instant. The true timescale will be less. 
For example, if fractures within the ejecta anneal, the rel- 
evant timescale is conductive cooling of a half-space (the 
ejecta layer) by an isothermal boundary condition (the well- 
mixed lake) until the heat flow into the bottom of the lake 
is less than heat loss at the top of the lake [ Turcotte & Schu- 
bert, 2002]: 

/ k( Tb -T) y i 

\EL vap + aT 4 J (ttk) V ; 

In this case, the conductive heat flow can only balance 
evaporative plus radiative losses for ~3 days for thermal 
diffusivity k = 10~ 6 m 2 s _1 : after this, an ice cover must 
form. 

Therefore, we are interested in spatially restricted (10° 
- 10 3 km) water sources which cease to emit vapor in 
timescales < 1 year. This is the domain of mesoscale 
modelling. We use the Mars Regional Atmospheric Mod- 
eling System (MRAMS), also used for entry, descent and 
landing simulations for the Mars Exploration Rovers, Mars 
Phoenix, and Mars Science Laboratory (Appendix A; Rafkin 
et al. [2001]; Michaels & Rafkm [2008]). MRAMS explic- 
itly resolves the size spectrum of dust and water ice aerosol 
for both cloud microphysics and radiative transfer, so it 
is well-suited for our cloud-forming numerical experiments 
[Michaels & Rafkin, 2008]. 

2. Orders of magnitude 

Order-of-magnitude reasoning suggests the possibility of 
precipitation close to Martian vapor sources. 

Consider liquid water at the surface of Mars: the atmo- 
spheric temperature is similar to today. The surface is not 
entirely frozen because drainage or filling stirs the lake and 
mechanically disrupts the ice cover, convection mines heat 
from an underlying hot layer to balance evaporative cool- 
ing, subsurface discharge outpaces evaporation, or fumaroles 
and gas-charged fountains inject vapor and small droplets 
directly to the atmosphere. The injection rate is approxi- 
mately [Emanuel, 1994] 

Q h ~C D \V a \{r* s -r b ) (3) 

where Co ~ 10~ 3 is a surface exchange coefficient 
(Emanuel [1994], p. 484), V a ~ 10 ms _1 is anemometer- 
level wind speed, r* ~ 0.5 is the near-surface vapor mix- 
ing ratio, and r b ~ is the background water vapor mix- 
ing ratio. To convert this to the vapor mixing ratio in air 
that enters the buoyant plume we require a length scale 



(a vertical distance over which the vapor is mixed) and 
a timescale (during which the vapor is injected). A rea- 
sonable length scale is the thickness of the subcloud layer 
AZ subc i ou d ~ 2 km. A reasonable time scale is the fetch 
timescale tf etch = Di ake /\V a [ ~ 2 hours for a 65km lake 
(where Di a ke is lake diameter). This gives a vapor mixing 
ratio in air that enters the buoyant plume, r ~ 0.01 = 6 Pa. 

Convection initiation is made more likely by Mars' low at- 
mospheric temperatures. Air containing 6 Pa vapor will be 
supersaturated with respect to ice when T < 227K [Hardy, 
1998]. Mars today has an atmospheric surface layer 20K on 
average colder than its ground surface temperature (Mars 
Climate Database v4.3, described by Lewis et al. [1999]; 
Millour et al. [2008]; henceforth Mars Climate Database). 
For example, at L s = 255° (near perihelion) the latitudinal 
maximum zonal mean daily maximum surface temperature 
is 304K, but the zonal mean daily maximum atmospheric 
surface layer temperature at this latitude is only 265K. This 
offset is due to the low atmospheric column density, be- 
cause this reduces radiative and mechanical coupling be- 
tween the atmosphere and surface. In the current climate, 
diurnal mean equatorial atmospheric temperature is below 
227K at all altitudes and all times of the year (Mars Climate 
Database). Because of these low temperatures, supersatu- 
ration of the vapor (and nucleation and growth via deposi- 
tion) occurs close to the ground. On Earth, the strongest 
lake storms are associated with rare cold-air outbreaks and 
large air-lake temperature contrasts [Markowksi & Richard- 
son, 2010]. On Mars, these conditions occur everyday. 

Mars' low atmospheric pressure promotes deep convec- 
tion (Figure 1). Because of the low volumetric heat capac- 
ity, we assume that the plume will accelerate upward until 
90% of the vapor has crystallized. Assuming that no precip- 
itation occurs, with a measured Mars adiabatic lapse rate 
of Tuiars ~ 1.5 K/km (Mars Climate Database) and the 
Clapeyron slope for ice at 230K, this occurs 12 km above 
cloudbase. 

During this ascent, the plume will have gained Convec- 
tively Available Potential Energy (CAPE): 

CAPE « PairQMars / — - (1 + n) dz (4) 

J C B 

where MC is the elevation of almost-complete condensa- 
tion, CB is cloud base, T is temperature within the plume, 
T" the environmental temperature, and n the ice mixing ra- 
tio [Rogers & Yau, 1989]. This assumes that there is no 
precipitation of ice out of the parcel. Approximating the ice 
crystallization as linear from at CB to complete at MC, 
we obtain T-T' = rL v /c C 02 ~ 30 K so T/T' ~ 1.15 upon 
complete crystallization, p, ~ 0.01 can then be set aside as 
negligible. Thus CAPE gained during ascent ~ 3000 J/kg 
and peak vertical velocity W ma x = \JlC APE / pair ~ 80 
m/s. We have ignored differential pressure gradient acceler- 
ation and compensating downward motions [Rogers & Yau, 
1989] in obtaining this result. The plume will continue to as- 
cend well above MC, but will decelerate as it entrains more 
ambient air and spreads to form an anvil cloud. The ascent 
timescale t ascend is 12 km / 0.5W ma x = 300s. Attention now 
shifts to the growing ice crystals. 

Although Mars air is cold enough at these altitudes to 
allow homogenous nucleation, we assume heterogenous nu- 
cleation occurs on dust. Equatorial dust opacity Td is typi- 
cally 0.01-0.05 during northern summer and 0.1-0.8 during 
southern summer [Liu et al, 2003]. The present day typ- 
ical low-latitude effective radius of dust is ~1.6/im [Wolff 
& Clancy, 2003]. A uniform dust density in the lower 30 
km implies a number density of 4 x 10 4 m~ 3 (4 x 10 6 m~ 3 ) 
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ice nuclei (considering geometric cross-section only and ne- 
glecting self-shadowing) for r = 0.01. At 20 km elevation (3 
g/m air; 0.03 g/m H2O), this yields an 'seeder' crystal ra- 
dius r of 58 ftm (or 13 /im for r = 1) - a minimum, in that it 
assumes all ice nuclei are consumed. Since crystallization is 
implicit in the updraft-velocity calculation, we do not assign 
a separate timescale to nucleation and early growth. 

How fast will crystals grow to precipitable size? Assum- 
ing that the 'seeder' crystals sink relative to the updraft 
and anvil cloud while scavenging both vapor and smaller 
droplets, the growth rate is just [Rogers & Yau, 1989] 



dR 

df 7 



EM 



AW 



(5) 



where E « 1 is collection efficiency, M ~ 0.03 g/m 3 H 2 
is cloud total water content at 20 km, pi = 910 kg/m 3 is 
ice density, and AW = 10 ms _1 is a characteristic sink rate 
relative to the surrounding vapor-laden air. This gives a 
growth time t growt h ~ 2000 s to precipitable size (assumed 
200 pm). In Earth thunderstorms this also requires tens of 
minutes. 

The crystals now begin to fall. At the low temperatures 
encountered at high altitudes on Mars, we expect the water 
to form hexagonal cylinders [Wallace & Hobbs, 2006]. The 
drag coefficient Cd for cylinders is ~1 over a wide range of 
Reynolds number Re, 10 2 - 10 5 [Tritton, 1988]. Therefore 
we obtain a terminal velocity of [Tritton, 1988] 



2g M arsPi^r 

0-5p a 



(6) 



which with r = 200 /im gives 40 ms" . Using the dy- 
namic viscosity of CO2 at 233K, p a = 1.2 x 10~ 7 Pa s, wc 
obtain Re — paU^r/pa ~ 200, sustaining the assumption of 
10 2 < Re < 10 5 . Fall time from anvil cloud height is then 
tfaii = Hpiume/uov ~ 500 seconds. 

The total lifetime of the vapor from release to precipita- 
tion as snow is t = t a scend+tgrow+tf a ll ~ 3000 s. During the 
entire process of ascent through the plume, crystal growth, 
and snow fall, the parcel has been blown sideways by the 
regional winds. Taking 20 m/s as a representative shear ve- 
locity at cloud-forming altitudes of ~ 20 km (Mars Climate 
Database), we find that snowfall will be roughly ~ 60 km 
downwind of source. This is small compared to the sizes of 
many Martian geographic features (craters, canyons, volca- 
noes). Therefore, provided latent heating powers a strong 
buoyant plume immediately downwind of the lake that lofts 
the released vapor to a height at which it will condense, we 
hypothesize that localized precipitation can occur on Mars. 

In order to test this hypothesis, we carry out a set of 
numerical experiments. 

3. Model setup 

Our numerical experiments used four nested grids with 
160 km resolution on the outermost (hemispheric) grid, in- 
creasing to 5.9 km on the innermost grid. The lake is cen- 
tered at 6.5S, 299E (the location of Juventae Chasma; paper 
2). 

For these simulations we prescribed flat land at 0m on all 
grids, with uniform albedo = 0.16, uniform thermal inertia 
= 290 kieffers [Mellon et al, 2000], and uniform roughness z 
= 0.03m (the value determined for the Mars Pathfinder land- 
ing site; Sullivan et al. [2000]). We introduced an isother- 
mal lake with albedo = 0.05, and constant roughness z = 
0.01m, which is close to the time-averaged roughness value 



obtained with the sea surface roughness parameterization 
of Eq. 7.21 of Pielke [2002]. Lake surface temperature is 
pinned to 278. 15K, with saturation vapor pressure accord- 
ing to Hardy [1998]. Constant lake surface temperature is a 
reasonable approximation if either (1) the lake is deep and 
well-mixed (e.g. for the 4km-deep lake at Juventae Chasma, 
cooling rate ~0.01 K day -1 for an evaporation rate of 2 
mm/hr if the lake is well-mixed), or (2) the lake surface 
is constantly refreshed by discharge of warm, perhaps gas- 
charged water from an aquifer [Harrison & Grimm, 2008; 
Bargery & Wilson, 2010]. In either case, the lake temper- 
ature will change more slowly than the (strongly diurnal) 
atmospheric response. The focus of this first study is to 
identify the steady-state response of the atmosphere to the 
lake perturbation, so a time-dependent treatment of coupled 
lake thermodynamics is not appropriate. We use the NASA 
Ames MGCM [Haberle et al, 1993] to provide atmospheric 
boundary conditions. To prevent extremely high water sub- 
stance mixing ratios and to allow metastable surface liquid 
water at +2 km above datum (the elevation of the Juventae 
plateau streams in Paper 2), we double atmospheric pres- 
sure relative to today in our simulation by doubling initial 
and boundary pressures supplied by the GCM. Runs are 
at L s « 270° (southern summer solstice). Runs were for 7 
days, except for the ref simulation which was extended to 
12 days to test for precipitation variability. Many of these 
parameters were varied in sensitivity tests: see Table 1 for 
a list of runs and parameters varied. More details on the 
model setup are provided in the Appendix. 

We carried out a dry run forced by these boundary con- 
ditions using flat topography but no lake. In the dry run, 
surface pressure is 1190-1270 Pa, 40% more than the satura- 
tion vapor pressure at 278. 15K. Surface temperatures range 
from 209-281K. Surface winds are from the NNE, backing to 
the N during the passage of each afternoon's thermal tide. 
Wind direction rotates anticlockwise with increasing height 
until, near cloudtop elevation (~ 30 km in the reference run), 
wind direction stabilizes at E to ESE (the subtropical jet of 
Mars' single Hadley cell). There are no significant day-to- 
day variations in the wind field. The mean 0-6km shear is 
22 m s -1 . 



4. Analysis of reference simulation 

For our reference simulation, referred to here as ref, wc 
introduce a square lake with sides ~ 65km. This is similar 
to the diameter of Mojave Crater, a young crater with a 
fluvially-modified rim [Mcewen et al, 2007; Williams & Ma- 
lin, 2008]). We begin the analysis by describing the time- 
averaged response. 

4.1. Time-averaged response 

Because the Mars daily average temperature is below 
freezing, the lake transfers sensible heat to the atmosphere 
on average. It also injects radiatively-active water vapor 
into the atmosphere, whose condensation at altitude re- 
leases latent heat and produces strongly-scattering water ice 
aerosols. Time-averaged surface temperatures up to 100 km 
downwind of the lake are raised by 8K as a result of these 
effects (Figure 2a). 

This number is small because the time averaged longwave 
heating (vapor greenhouse and cloud reradiation) and short- 
wave cooling (ice scattering) due to the lake nearly cancel. 
Time-averaged longwave forcing peaks at +106 W/m 2 im- 
mediately downwind of the lake, which is where vapor col- 
umn abundance is greatest. Time-averaged shortwave forc- 
ing peaks at -83 W/m 2 about 20km further SSE, which is 
where ice column abundance is greatest. Net radiative forc- 
ing is ~ +30 W/m 2 over the lake, but typically between -3 



X - 4 



KITE ET AL.: LOCALIZED PRECIPITATION ON MARS 



and +15 W/m in the area of greatest ice precipitation (Fig- 
ure 2b). Because of the near-uniform northerly wind, we can 
understand the spatial structure of the net radiative forcing 
in terms of the timescales needed for vapor to ascend and 
precipitate out:- net forcing is positive just downwind of the 
lake, where vapor is still ascending to condensation level; it 
is modestly positive or slightly negative while vapor is be- 
ing converted into scatterers; and it becomes positive again 
when those scatterers have had time (equivalently, distance 
South) to precipitate out. 

A vertical slice through the atmosphere shows three com- 
ponents to the time-averaged atmospheric temperature re- 
sponse. Atmospheric temperatures in a thin boundary layer 
above the lake are warmed by 30K. The planetary bound- 
ary layer is thickened downwind of the lake because of the 
increased turbulence associated with the lake, so tempera- 
tures rise by 4K in the part of the atmosphere that is in- 
cluded in the boundary layer as a result of the presence 
of the lake. Most importantly, a narrow plume of 5-10K 
increased temperatures extends 10-15 km above the lake. 
This thermal plume corresponds to the latent heat released 
by the lake- induced storm. Atmospheric temperatures never 
exceed 273K, so liquid water droplets are never stable. Su- 
percooled water droplets are not included in our simulations, 
and could only form very close to the lake surface. 

Strong low-level convergence results from the release of 
latent heat (Figure 2d). Because the lake effect creates its 
own convergence, this suggests that precipitation should be 
insensitive to changes in the regional windfield. 

The total water column abundance (thin lines in Figure 
2e) shows the narrow extent of the weather system induced 
by the lake: its core is similar in horizontal extent to a 
terrestrial thunderstorm. The contours of (ice/total water) 
fraction (thick lines in Figure 2e) are extended to the west 
of the lake because the ice-rich uppermost levels of the cloud 
are affected by the Easterly subtropical jet. Resublimation 
within this jet lowers the ice fraction with increasing dis- 
tance to the west. 

4.2. Spatially-averaged time dependencies 

The area of peak water-ice precipitation is immediately S 
of the lake. Within a square with sides ~65km immediately 
S of the lake, the daily temperature cycle is regular, with 
little variability. Water vapor and cloud blanketing raises 
nighttime surface temperature by up to 18K relative to the 
dry run, but the net increase in daytime surface temper- 
ature due to the lake is small or negative because of ice- 
particle scattering (Figure 3a). Despite ice-particle scat- 
tering, afternoon surface temperature exceeds 273K in the 
area of peak water-ice precipitation, so snow falling onto 
bare ground during the afternoon will melt. Snow falling 
onto ground that is cooled by the increased albedo of snow 
that has fallen during the night may or may not melt, de- 
pending on the grainsize, thickness, and dust content of the 
nighttime snow layer [Clow, 1987]. 

To track water and ice mass budgets, we average over 
a square 400 km on a side which contains the lake but is 
centered 100km S of the lake in order to enclose the cloud. 
The time dependence of the water substance mass budget is 
dominated by a strong afternoon peak in atmospheric wa- 
ter vapor (~2.5 x the predawn minimum of ~9 Mton) at 
a time when atmospheric temperatures are highest (Figure 
3b). The mass of water ice in the atmosphere is indepen- 
dent of time-of-day and averages 10 Mton during sols 5-7. 
The majority of the water vapor injected into the system 
precipitates as water ice during the night (see §4.4). Again, 
notice that atmospheric temperature does not exceed 273K, 
so liquid water aerosol is never stable. 

4.3. Structure of the buoyant plume 

Latent heat release drives plume ascent. Figure 4 is a 
predawn snapshot: water ice mass ratio in the plume core 
exceeds 1%. At 17.5 km elevation, more than 95% of water 



substance is in the condensed phase and vertical velocities 
reach 54 m s _ . There is now little energy to be gained 
from further condensation, so the plume slows and broad- 
ens. Sublimation, entrainment, and especially precipitation, 
all lower the ratio of ice to vapor. The plume overshoots its 
equilibrium level, peaking near 35km. 

At a given altitude, plume updraft velocities are strongest 
just before dawn and weakest in early afternoon. This di- 
urnal cycle in plume behavior corresponds to three related 
changes in the state of the atmosphere just upwind of the 
lake. Firstly, the Planetary Boundary Layer (PBL) pinches 
and swells during the diurnal cycle. The greater depth and 
intensity of turbulent mixing during the day leads to more 
entrainment of moist air by the ambient atmosphere, where 
it does not contribute to the plume. Higher air temperatures 
disfavor crystallization at low altitudes during the day. To- 
gether, these two changes allow more vapor to escape during 
the day, by advection downwind within the thickened PBL, 
and not contribute to the plume. Finally, during the day, the 
excess of lake temperature over land temperature is small, 
and convergence is weak, so the shear velocity tt* is small. 
Therefore, relatively little vapor mixes from the lake surface 
boundary layer into the atmosphere. During the night, the 
greater land-lake temperature contrast is associated with 
stronger convergence. More vapor is mixed above the sur- 
face layer and entrained by the plume. Allowing for dynamic 
lake surface roughness (zo ~ u,) increases the sensitivity of 
plume velocity to time-of-day, which confirms the effect of 
surface-layer dynamics on the diurnal cycle. 

4.4. Precipitation 

Precipitation is strongly peaked just downwind of the 
lake. Figure 5a shows that time-averaged snowfall is ev- 
erywhere < 0.006 x its peak value at distances > 150 km 
from that spatial peak. Peak snowfall (Figure 5b) shows a 
similar, but more noisy, pattern, with peak precipitation > 
4.5 mm/hr only in a small area beneath the buoyant plume. 
Precipitation rates increase by a factor of 4 during the night, 
with a rapid decline during the morning to a lower, stable 
afternoon rate (Figure 5c). 

Diurnally-averaged precipitation is steady in location over 
the 12 days of our extended ref simulation. The total pre- 
cipitation rate increases slightly with time. 

4.5. Precipitation efficiency 

Figure 6 and Table 2 show the fate of released vapor at the 
end of sol 6 in ref. The majority of the vapor precipitates 
<200 km from the lake. Of the vapor that reaches distances 
>200 km from the lake, the majority is in the form of vapor 
at the end of the simulation. 

After 6 days of the ref run, the outermost (hemispheric) 
grid contains ~300 Mt more atmospheric water than the dry 
run (please refer to Table 1 for details of runs and param- 
eters). This is equivalent to a global surface liquid water 
layer 2pm thick, or a lake depth of 7 cm. This is radia- 
tively unimportant on the global scale, and less than the 
current global average (~17 pr /im in the northern hemi- 
sphere, ~9.5 pr /im in the South: Smith [2002]). The added 
vapor will presumably precipitate as ice on the winter pole, 
which is outside the space and time limits of our simulation. 
If all vapor precipitates in a single season on a polar cap of 
area 10 km 2 , it would form a layer 0.3 m thick. Layers of 
this thickness could be resolved by the HiRISE camera on a 
gently-sloping exposure. 

Figure 6 understates the fraction of water that precipi- 
tates locally if the conditions that maintain a liquid surface 
are maintained indefinitely. This is because some of the wa- 
ter that is in the atmosphere at the end of sol 6 will precipi- 
tate locally. In an extended run (ref (sol 12) in Figure 6), 
the fraction of water that precipitates locally is increased. 



KITE ET AL.: LOCALIZED PRECIPITATION ON MARS 



X - 5 



5. Sensitivity tests 

5.1. Vertical resolution 

We carried out a sensitivity test to determine model va- 
por release as a function of vertical resolution. The penalty 
of increasing vertical resolution scales as iVlog(iV), so runs 
were for only \ sol starting at ~ 12 noon Local Solar Time. 
Vapor release rate is moderately sensitive to model vertical 
resolution. Increasing surface layer thickness by a factor of 
30 decreases vapor release by 35% (Figure 7). Therefore, our 
modeled afternoon vapor release (and precipitation) rates 
are probably underestimates. (In an earlier report (Kite et 
al. [2010]), we stated the results of a sensitivity test carried 
out at night, for which the sensitivity has similar magnitude 
but opposite sign.) 

5.2. Horizontal resolution 

We carried out a simulation, hires, that added an inner 
grid with 2.0 km horizontal resolution, (hires was initial- 
ized from ref output after 2 days) . Overall storm structure 
was similar, although secondary plumes developed in addi- 
tion to the main plume. Peak vertical velocity increased 
by 34%, peak time-averaged vertical velocity by 57 %, and 
cloud height by 12%, relative to ref. However, when aver- 
aged to the lower resolution of ref , these differences decrease 
to 16%, 26%, and 12% (unchanged), respectively (Figure 
8). We conclude that our main results are not affected by 
horizontal resolution. Simulations which use resolutions in- 
termediate between eddy-resolving and the mesoscale are 
known to suffer from artifacts caused by aliasing of barely- 
resolved eddies by the grid cell size. Given that we observed 
2 Ax noise in hires, we consider our reference run to be more 
accurate. 

5.3. Size of idealized lake 

Small lakes in our simulations are unable to drive deep 
convection, and have weaker localized precipitation. We 
modeled square lakes with areas of ~35 km 2 (s), ~300 km 2 
(m), ~1700 km 2 (l), and ~29000 km 2 (xxl), in addition to 
the ~4000 km 2 ref simulation. 

Deep moist convection was not observed in the s and m 
simulations. To verify that this was not an artifact of in- 
sufficient model resolution, we ran a 2km-resolution nested 
grid on m. Although peak vertical velocities did increase in 
the nested grid relative to the default-resolution model, lake- 
sourccd vapor did not reach altitudes much greater than the 
planetary boundary layer and, as in the default-resolution 
model, was passively advected downwind. Given our par- 
ticular choice of atmospheric boundary conditions (present 
day orbital conditions, low latitude, and southern Summer), 
lake size > 10 3 km 2 is required for deep moist convection. 
If vapor is funnelled to a single central plume of constant 
radius and entirely condensed, we would expect plume ver- 
tical velocity to scale linearly with lake area. Figure 9 shows 
that this expectation is borne out for lakes of size 10 2 ' 5 km 2 - 
10 3 ' 5 km 2 , but overestimates lake storm convective intensity 
for the largest simulation. This xxl simulation is unusual 
because it has a much broader area of strong upwclling than 
the smaller lakes, which explains why its convective veloci- 
ties are not as high as expected. In the smallest lake mod- 
eled, the lake-associated updrafts are so weak that they were 
difficult to seperate from everyday planetary boundary layer 
convection. 

Because lake-induced convergence (Figure 2d) efficiently 
funnels lake-released vapor into one buoyant plume, plume 
velocity and plume height increase with lake size (Figure 8). 
Greater plume heights promote ice formation, and the ice 
fraction of atmospheric water increases with lake size (colors 
in Figure 6). This is a self-sustaining feedback, because ice 
formation provides energy for plume ascent, which in turn 
creates low pressure above the lake and drives convergence. 



Larger lakes inject proportionately less water into the 
global atmosphere. More than | of the water released in 
the s,m and 1 simulations is in the atmosphere at the end 
of the simulation. <14% is in the form of ice. In the xxl 
simulation, the fraction of released water that is precipi- 
tated locally is >90%, and most of the atmospheric water 
is ice aerosol that would probably precipitate locally if the 
lake surface froze over. The line source has a similar lake 
area to the square xxl run, and behaves similarly in having 
a small atmospheric water fraction and a large ice aerosol 
fraction. This suggests that this trend of increased local- 
ized precipitation fraction with increased lake area does not 
depend on lake geometry. 

5.4. Surface roughness parameterization 

Our default run assumes fixed lake surface roughness, but 
water surface roughness increases with wind speed. A more 
realistic approximation is that of Garrett [1992] as given by 
Eq. 7.21 of Pielke [2002]: 

z = 1 x 10~ 4 + (0.01625u 2 /p) (7) 

where we have added the first term to maintain numerical 
stability. This remains an adequate fit to the much larger 
datasets now available [Zeng et al., 1998]. The time-average 
lake roughness with this parameterization is 0.0126 m, and 
is a factor of three higher during the night when near-surface 
winds are stronger. 

Our time-averaged results are insensitive to this more 
accurate lake surface roughness parameterization (Table 2; 
Figure 6), although the very strong plumes seen during the 
night (when vapor injection rate is highest) increase both 
plume height and updraft velocity (Figure 8). 

5.5. Line source 

Our line source is intended to sketch an outflow channel 
during a channel-forming flood. It is a N-S oriented, straight 
trough of depth 1.5 km, floor width 30 km, wall slope 0.13, 
and length 880 km. The floor is flooded. Of the resulting 
snow, 23% falls back into the lake and a further 70% falls 
within 100 km of the edge of the trough. The diurnal cycle 
consists of a strong, steady, spatially continuous line storm 
on the W edge of the trough during the night, and a clumpy, 
broken line of weak updrafts some distance E of the trough 
during the day. Time-averaged vertical velocity and precip- 
itation fields do not show this dumpiness: intead, there is a 
trend of monotonically increasing vertical velocity and pre- 
cipitation rates toward the north, because the background 
wind field advects the upper parts of the cloud toward the 
north. 

We do not consider the influence of the flowing water on 
surface roughness, nor the drag of the flowing water on the 
atmosphere. 

5.6. Season 

In southern winter (also aphelion season) the sign of the 
Hadley circulation reverses. In a run at Ls = 90°, the GCM 
boundary conditions produced ESE-directed time-averaged 
winds at altitudes below ~ 15 km, and SW-directed time- 
averaged winds at higher altitudes. The highest snowfall 
rates were just ESE of the lake, reflecting steering of the 
plume by regional winds. Total localized precipitation was 
the same as in ref to within 7%. 

5.7. Latitude 

An cyclonic circulation does not develop in our ref run. 
To determine if Coriolis effects can produce an cyclone at 
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higher latitudes, we ran a test at 45° S latitude. As expected 
given the small size of the lake, the lake-driven circulation 
is too weak to restructure the background wind field and a 
lake-induced standing cyclone does not form. 

The total localized precipitation in the midlatitude run is 
only | of the total localized precipitation in the equatorial 
runs. Since vapor release is similar, vapor is being con- 
verted to localized precipitation with a smaller efficiency. 
Peak time-averaged column ice fractions are < 40% in the 
midlatitude run, compared to > 70% in the reference run 
(Figure 2e). This is because both runs are close to southern 
summer solstice. The southern midlatitude site is in sun- 
light for 2/3 of each day, and has air temperatures higher 
by 20K on average than the equatorial site. Because of the 
lower supersaturations, water in the midlatitude plume must 
be lifted 5 km higher to obtain a given ice fraction than in 
the equatorial plume. Therefore, more vapor escapes to the 
regional atmosphere. 

5.8. Faint young sun 

Many channel networks on Mars formed when the Sun 
was fainter. The solar luminosity L at the Hesperian- 
Amazonian boundary was 0.78-0.85 x present-day solar lu- 
minosity Lnow, and L at the Noachian- Hesperian bound- 
ary was 0.75-0.77 x Lnow [Bahcall et al, 2001; Hartmann, 
2005] . The uncertainty is due to differences among the mod- 
els that map crater density onto absolute age: the solar evo- 
lution model has much smaller error bars. The Juventae 
plateau channel networks (Paper 2) could be as old as Hes- 
perian [Le Deit et al, 2010], so we ran a test at 0.75 Lnow- 
As well as being less luminous, the young Sun was also 1.5 
% cooler [Bahcall et al, 2001]. We ignore the resulting small 
shift in the solar spectrum and simply reduce the flux at all 
wavelengths by 25 %. 

Colder air temperatures under the faint young sun lower 
the cloud base by ~ 5 km. Peak updraft velocities are >50 
m.s. in both simulations, but occur at lower elevations in 
the faint young sun simulation. Because supersaturations 
are higher at all colder altitudes, ice growth is favored and 
the fraction of atmospheric water that is ice increases from 
32 % in ref to 47 % in lo_sun. 

We conclude that localized precipitation is favored by re- 
duced solar luminosity. However, this is not true for runoff 
generation and channel formation (§6). 

5.9. Paleoatmospheric pressure 

Figure 1 suggests that localized precipitation will not oc- 
cur if atmospheric pressure is greatly increased. To test 
this, we increased pressure to lOx Present Atmospheric Level 
(PAL) on all the mesoscale grids (run hipress). Deep 
moist convection is suppressed, uplift velocities are reduced, 
and comparatively little vapor reaches cloud-forming alti- 
tudes. No lake-sourced ice is found above 10km, and little 
ice aerosol forms (1.6 % of atmospheric water mass on the 
inmost grid, versus 32 % in ref). 

Precipitation is reduced by 20 %, and preliminary runs 
at 20x PAL lead us to expect that further increases in pres- 
sure will greatly reduce localized precipitation. However, the 
mismatch between the forces driving winds on the GCM (at 
lx PAL) and and the forces driving winds on the mesoscale 
grids at lOx PAL is already severe. Higher-pressure simula- 
tions will benefit from higher-pressure GCM runs. 



6. A probabilistic model 
melting on ancient Mars 



of snowpack 



For a lake lifetime of 1 Earth year (similar to expected 
chaos outflow durations: Andrews- Hanna & Phillips [2007]; 
Harrison & Grimm [2008]), and assuming that seasonal ef- 
fects on the storm are minor, ~2 m of snow is predicted to 



accumulate downwind of the lake by our reference simula- 
tion, assuming a snowpack density of 350 kg/m 3 . Will it 
melt? 

We expect the lifetime of snowpack against sublimation 
and wind erosion to be > O(l) year, so that the annual- 
maximum temperature is that which is relevant for snow 
melting. For predicted snowpack depths 0(l)m, this resi- 
dence time is supported by GCM simulations (e.g., Table 1 
in Madeleine et al. [2009]). A simple energy argument gives 
the same result - 



tsubli 



d snow p snow L s 
FsunT/fd 



(8) 



which with snow depth d sn0 w = 1 m, snow density p sn ow 
= 350 kg/m 3 , solar flux Fs un = 300 W/m 2 , fraction of solar 
energy used for sublimation r\ — 0.1, day fraction fd = 0.5, 
gives tsublime = 10 Mars year. With surface pressures on 
the plateau of ~1000 Pa, maximum modelled surface shear 
stress is not sufficient to initiate motion of snow that is rest- 
ing on the ground [Greeley & Iversen, 1989], so snow cannot 
be blown away. 

Atmospheric temperatures in our model are too low for 
rain to be stable, so precipitation falls as snow. This is 
true even for runs (not shown) in which initial surface tem- 
perature is that of impact ejecta: for a low-pressure, ra- 
diatively thin atmosphere similar to today, energy transfer 
from the surface to the atmosphere is slow and inefficient. 
If localized precipitation in a thin atmosphere similar to to- 
day can account for channel formation, it must be shown 
that surface liquid water runoff can occur on <5° slopes 
in a thin atmosphere similar to today. Midlatitude chan- 
nels at Acheron ~0.08 Ga and Lyot < 1.5 Ga formed on 
<5° slopes sufficiently recently that it's unlikely that the 
atmosphere when they formed was much thicker than it is 
today [Dickson et al, 2009; Fassett et al, 2010]. This of- 
fers empirical support for the snowpack runoff hypothesis. 
However, Amazonian midlatitude channels are not common, 
and formed in unusual circumstances (such as nearby steep 
slopes, and glaciers) [Fassett et al, 2010] that may not hold 
for low-latitude, localized channels. In addition, some of the 
observed low-latitude channels (e.g, those described in Pa- 
per 2) formed > 3.0 Ga when solar luminosity was <80% 
that of today's Sun [Bahcall et al, 2001]. A model of the 
likelihood of snowpack melting is needed. 

In general melting probability will depend on orbital ele- 
ments including obliquity and precession, latitude, age (via 
solar luminosity), material properties, and the lifetime of 
snowpack against sublimation losses to planetary cold traps. 
For age >20 Mya, deterministic chaos makes orbital ele- 
ments unreliable: as our interest is in ancient Mars our 
model is therefore probabilistic. We model temperatures 
at the equator - 64% of the sedimentary rocks on Mars are 
within 10° of the equator (Malin et al. [2010], figurel6.txt 
in their online supporting data), and our modeling target in 
Paper 2 is at 5° S. 

a. A a simple energy balance. The energy balance of a 
snow-bearing surface on Mars is [Dundas & Byrne, 2010] 



dU 
dt 



(1 - a)SWi +eLWi -eaT 4 - SH - C net - L s 



duiw 

'~dT 

(9) 



where the left-hand side is gain of energy by the surface 
layer, a is surface albedo, SW \, is sunlight, e ~ 0.98 for 
ice is thermal emissivity, LWi is the greenhouse effect, a 
is the Stefan-Boltzmann constant, T is surface temperature, 
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SH is net sensible heat flux from the surface to the atmo- 
sphere, Cnet is conductive heat flux from the surface into 
the ground, L suo is the latent heat of sublimation of water, 
and d" t m is the sublimation rate. 

The minimum sublimation rate (in the absence of wind) 
is given by [Hecht, 2002; Dundas & Byrne, 2010] 



am ™ n 1 a A n M A P a " 
—— = 0.14Ar)p ave D 

Ot \ \ Pair 



QMars 



\ D a i r / 



1/3 



(10) 



where An is the difference between atmospheric and sur- 
face gas water mass fractions, p ave is the atmospheric den- 
sity, D air is mass diffusivity, Ap air /pair is the normal- 
ized density difference between the moist near-surface at- 
mosphere and ambient atmosphere, and v air is kinematic 
viscosity. We parameterize the temperature-dependency of 
Vair and D a i r following Dundas & Byrne [2010]. 

Melting occurs if ^ > at T = 273. 15K. We set up 
optimistic conditions for melting, and ask: with what prob- 
ability do orbital conditions permit melting? 

Evaporative cooling is handled following [Dundas & 
Byrne, 2010]: we obtain 165W assuming 70% humidity, at- 
mospheric pressure 1220 Pa, and zero wind. The more so- 
phisticated surface-layer theory of Clow & Haberle [1990] 
gives an evaporative cooling rate >2 times less than in Dun- 
das & Byrne [2010] (and >4 times less than in Ingersoll 
[1970]) for roughnesses «o <0.3 mm that are appropriate for 
polar snow [Brock et al., 2006]. Therefore, our parameteri- 
zation may overestimate evaporative cooling by free convec- 
tion. We assume still conditions, so that cooling by wind- 
driven turbulent fluxes (forced convection) is negligible: this 
assumption is favorable for melting. We also neglect the 
turbulent flux of sensible heat SH, which is a factor of >10 
smaller than evaporative cooling for these conditions [Hecht, 
2002]. The noontime conductive heat flow C„ e t can be pa- 
rameterized assuming the temperature falls to the diurnal 
average temperature at a depth equal to twice the diurnal 
skin depth. Assuming a diurnal-average temperature during 
the hottest season of 220K, this gives 



80°), eccentricity (e, 0.0 — > 0.175), longitudes of perihe- 
lion (L p ) and solar longitudes that are sampled by Mars 
[Laskar et al, 2004]. We use a ID thermal model to calcu- 
late temperature T within snowpack at a range of seasons, 
ignoring the greenhouse effect LW\.. T is solved for by ma- 
trix inversion for the linear (conduction) part and iteratively 
for the nonlinear contribution of radiation and evaporative 
cooling in the topmost layer. The timestep is 12 seconds. 
Evaporative cooling is handled following [Dundas & Byrne, 
2010] with optimistic assumptions as in the simple energy 
balance. Each run is initialized with surface temperature 
T SU rf at instantaneous thermal equilibrium with incoming 
sunlight, decaying with an e-folding depth equal to the di- 
urnal skin depth d = U (kp) / (Pc v ty) to an energy-weighted 
time-averaged equilibrium temperature at depth. (Here, k 
is thermal conductivity, p density, P the length of 1 sol in 
seconds, and c p specific heat capacity of the subsurface.) We 
then integrate forwards in time for multiple sols, but with 
seasonal forcing held constant, until sol peak T sur f has con- 
verged. Models that neglect the greenhouse effect provide a 
good approximation to Mars' observed surface temperature. 

Atmospheric contributions LWl are parameterized as a 
time-independent greenhouse temperature increase, AT a t- 
The purpose of adding the atmospheric contribution in post- 
processing is to allow different models of the past atmo- 
sphere to be compared to each other without rerunning the 
ensemble of ~10 5 column models. Typical present-day val- 
ues are AT at = 5-10K [Read & Lewis, 2004]. This aver- 
ages over strong nighttime warming, which is irrelevant for 
melting, and weak afternoon warming, which is critical to 
melting. Setting AT a t to the time average of present-day 
greenhouse warming will therefore lead to an overestimate 
of melting. Instead, we want to know the value AT a t at 
the peak temperatures relevant to melting. We make the 
approximation 



AT a - : 



' surf, max 



surj,max 



LW Xt 



(12) 



_ ,9T _ , 273.15 - 220 

(-/net — _ — K 

dz 2y/nP/n 



(11) 



For thermal conductivity k = 0.125 W/m/K for snow 
[Carr & Head, 2003], sol length P = 88775.204 s, and ther- 
mal diffusivity n = 2.04 x 10~ 7 m 2 s _1 appropriate for low- 
density snow, we obtain C ne t = 40.3W/m 2 . 

The greenhouse effect LW 4- at the time of peak sur- 
face temperature is assumed to be 55W/m 2 (Mars Cli- 
mate Database). We do not consider the additional green- 
house forcing from the lake storm, nor from a possible past 
stronger greenhouse effect. 

For dusty snowpack with albedo 0.28, the resulting min- 
imum on SW is 640W. 

We calculate the maximum equatorial luminosity for 
all ages, seasons, and orbital conditions. The age- 
dependent solar luminosity is taken from the standard so- 
lar model of Bahcall et al. [2001]. We then weight the 
annual-maximum luminosity results by the age-dependent 
probability densities for obliquity and eccentricity from 
Laskar et al. [2004], assuming that obliquity and eccen- 
tricity are not strongly correlated. (In the 0.25 Gyr so- 
lutions provided by J. Laskar at http://www.imcce.fr/ 
Equipes/ASD/ insola/mars/mars.html, obliquity and ec- 
centricity are not strongly correlated.) 640W is exceeded 
with 21% probability at 1 Gya, but with only 0.5% proba- 
bility 3 Gya (Figure 10). 

b. A ID column model. Next we calculate temperatures 
within snowpack for the full range of obliquities ((f), 0° — > 



The Mars Climate Database Mars Year 24 simulation 
shows ATat 5-7K at low latitudes. This excludes the dust 
storm season, for which the neglect of atmospheric scatter- 
ing leads to an overestimate of the net atmospheric contri- 
bution. We therefore adopt 6K as our greenhouse forcing. 
Additional greenhouse forcing due to doubled CO2 is small 
and safely neglected [Wordsworth, 2010]. 

Material properties for snowpack are taken from Carr & 
Head [2003]: we assume a snowpack density of 350 kg/m 3 , 
thermal conductivity of 0.125 W/m/K, and specific heat 
capacity c p of 1751 J kg -1 , and neglect the temperature de- 
pendence of k and c p . We consider dirty H2O snow albedos 
between 0.28 and 0.4. The smaller value is the present-day 
albedo of Mars' dust continents [Mellon et al., 2000], and the 
larger value is the upper end of the envelope that best fits 
the observed seasonal dependence of ephemeral equatorial 
snow on present-day Mars [Vincendon et al, 2010]. (Paper 
2 provides a more detailed discussion of snowpack albedo). 

Maximum-temperature results are interpolated onto a 
finer mesh in orbital parameters (<j>, e, L p ). We again weight 
the results using the probability distributions of Laskar et 
al. [2004]. Results for albedo = 0.28 and latitude = 0° are 
shown in Figure 11. Temperatures greater than 273. 15K arc 
not realistic, because snowpack temperature will be buffered 
by the latent heat of melting. The agreement between the 
melting probabilities obtained from the ID column model 
and the melting probabilities obtained from the simple en- 
ergy balance argument is excellent (Figures 10 and 11). 
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Even with the weak, present-day greenhouse effect, temper- 
atures high enough for melting occur during 4% of the years 
(0.1% for albedo 0.35). With an additional 6K of green- 
house warming, melting occurs during 56% of years (12% 
for albedo 0.35). 

Discussion. The exact melting probabilities are unlikely 
to be correct, but we draw three lessons from these simple 
models: (1) The probability of low-latitude melting is cer- 
tainly higher for instantaneous emplacement of snow than 
for gradual (e.g, orbital equilibrium) emplacement of snow, 
and probably much higher. That is because gradual em- 
placement of annually-persistent snow at the equator is only 
possible at high obliquity (~50% of the time), while in- 
stantaneous emplacement can occur at any time. In ad- 
dition, gradually-emplaced snow will accumulate preferen- 
tially in areas where it is most stable: for example, shad- 
owed areas, adverse steep slopes, and especially chasma and 
crater interiors where shadowing reduces peak temperature. 
Craters at 10° S in Sinus Sabeus appear to contain mantled, 
atmospherically-emplaced ice deposits, providing direct ev- 
idence that equatorial ice on Mars does preferentially accu- 
mulate in crater-interior cold traps [Shean, 2010]. The most 
compelling evidence for low-latitude glaciation on Mars is on 
the flanks of the Tharsis Montes, which have the lowest at- 
mospheric pressures on the planet [Forget et al, 2006]. Liq- 
uid water exposed on the Tharsis Montes today would boil 
internally, leading to very rapid evaporative losses. Taken 
together, these observations suggest that localized precipi- 
tation is much more likely to melt than orbital-equilibrium 
precipitation. 

Today's orbital conditions would lead to melting of snow- 
pack on the equator (diamond in Figure 10). However, on 
today's Mars, snow is almost entirely restricted to high lati- 
tudes. Snow and ice on Mars do not melt in general because 
the area of transient liquid water stability changes on slow, 
orbital timescales O(10 4 ~ 5 ) yr. This allows time for ice to 
be removed (vertically or latitudinally) from the advancing 
zone of increased saturation vapor pressures, towards cold 
traps where melting cannot occur. Therefore the areas of 
surface ice and of transient liquid water stability rarely in- 
tersect. Theory, experiments, and the lack of evidence for 
transient liquid water at the Phoenix landing site confirm 
this [Schorghofer, 2010; Hudson et al., 2007; Mellon et al., 
2009]. Steep (>20°) gullied slopes within mid- and high- 
latitude craters are an exception [Morgan et al, 2010], as- 
suming that the gullies formed from liquid water flows. In 
contrast, localized precipitation delivers snow on timescales 
O(10 _1 ) yr to a location uncorrelated with transient liquid 
water stability, so a wider range of orbital conditions will 
then allow melting. 

In the case of impact-induced precipitation, snow falling 
on hot ejecta will melt regardless of orbital conditions. 

(2) Because of the wide range of possible orbital states, 
the probability of melting is not large, but neither is it zero. 
The distribution of melting events (T sur f >273K) with or- 
bital parameters is shown in Figure 12. Relative to the pdf of 
obliquity, melting events are more probable at low obliquity 
than high; they are much more probable when perihelion 
occurs at equinox than when perihelion occurs at solstice; 
and they occur exclusively at moderate-to-high eccentricity. 
These results may be understood as follows. Moderate ec- 
centricity is needed to offset the reduced solar luminosity. 
For example, at 3.0 Ga and e = 0.11, perihelion insolation 
is eqiuivalent to insolation on a circular orbit today. Peri- 
helion near equinox is needed to align the season when the 
noontime sun is highest at the equator with the peak solar 
flux. Increasing obliquity narrows the seasonal window dur- 
ing which the sun rises high enough in the equatorial sky 
for melting to occur. Therefore, the probability of aligning 
with perihelion is reduced and so the chance of a melting 
event at high obliquity is lower. 

This distribution is almost exactly a ranking in peak solar 
luminosity, and is not sensitive to our model assumptions. 



(3) Using slightly less optimistic assumptions in the en- 
ergy balance approach (humidity = 0.25 and albedo 0.35), 
the minimum on SW is 882W. This is unachievable at Mars. 
For comparison, the sunlight striking the top of the atmo- 
sphere above North Greenland, 80N, at summer solstice dur- 
ing the melt season is 730W. The strength of the evaporative 
cooling term at low pressures strongly suppresses snowpack 
melting. Evaporative cooling may not be relevant for un- 
derstanding the melting of small quantities of snow that is 
in contact with non-volatile soil, and small quantities may 
be all that is necessary to redistribute mobile elements and 
to form crusts in the global soil [Amundson et al., 2008; 
Arvidson et al, 2010]. However, this cannot account for 
runoff generation and channel formation. Adding CO2 pres- 
sure quickly damps evaporative cooling, but also increases 
the magnitude of sensible heat transfer to the atmosphere. 
There is a minimum in cooling at intermediate pressures 
[Hecht, 2002], but this may suppress localized precipitation. 

For Early Mars, this dilemma can probably only be re- 
solved by an early orbital state with a smaller semimajor 
axis or smaller eccentricity than the current one; a multi- 
bar CO2/H2O atmosphere; or adding non-C02 greenhouse 
gases. In summary, our models do suggest a requirement 
for a different global climate state if localized precipitation 
is to lead to substantial runoff on gentle slopes on Early 
Mars, but this is only necessary to bring temperatures back 
up to the level of contemporary Mars: melting probabili- 
ties with present-day solar luminosity are large (Figure 13). 
To counteract the effect of the Faint Young Sun on mean 
Mars surface temperature at 3.5 Gya, the required green- 
house warming is approximately 



AT 



fys 



(1 - 0.77 a25 )T M ar 



13K 



(13) 



(with TjMars ~ 210K; the luminosity at 3.5 Ga is from 
Bahcall et al. [2001]). This is a relatively modest require- 
ment. It does not hold for impact-induced precipitation, 
because snow falling on hot ejecta will melt provided that 
atmospheric pressure is above the triple point. 

7. Prospectus and discussion 

Our model shows lake storm ice forming at terrestrial 
cirrus cloud temperatures >200K, which because of the 
higher water vapor mixing ratio is warmer than modern 
Mars cloud forming temperatures <185K. Recent laboratory 
experiments [Iraci et at, 2010] show that ice nucleation at 
T < 185K requires greater supersaturations than previously 
thought, while confirming that the MRAMS-CARMA pa- 
rameterisation of the critical supersaturation at the higher 
ice nucleation temperatures observed in our simulation is 
adequate. 

We have focussed on ephemeral open-water lakes in this 
study, but this is not the only way of getting water into 
the atmosphere. Fumaroles and rootless cones inject vapor 
to the atmosphere, but their efficiency has not been as well 
studied as that of a wind-stirred lake surface [Zeng et al., 
1998]. Gas-charged fountains may inject vapor during the 
early stages of outflow channel formation [Bargery & Wil- 
son, 2010]. The mass flux we obtain in even the smallest 
simulation greatly exceed the greatest spring discharges on 
Earth. Leads in lake ice may allow deep convection to con- 
tinue after the majority of the lake surface has frozen over. 

7.1. Geomorphic paleobarometer? 

Cold lakes at 1 bar on Earth do not produce vigorous 
storms, with convection and precipitation in a dry back- 
ground atmosphere (Figure 14); cold lakes at 12 mbar on 
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Mars do, at least in our model. Therefore there is a thresh- 
old atmospheric pressure, between 12 mbar and 1 bar, below 
which cold lakes on Mars can induce convective storms and 
localized precipitation (orange line in Figure 1). We specu- 
late this pressure is ~10 2 mbar, because this is the pressure 
at which a 273K lake on Mars produces the same buoy- 
ancy flux as 299K on Earth. 299K is the lower limit of the 
tropical sea surface temperatures that on today's Earth are 
associated with deep convection [Emanuel, 1994]. A prelim- 
inary run at 20x present-day pressure (~ 120 mbar) showed 
significant suppression of the buoyant plume, with limited 
localized precipitation. If a past landform on Mars could 
be conclusively shown to result from localized, lake-induced 
precipitation, that would suggest a low atmospheric pressure 
(< 10 2 mbar) at the time it formed. We investigate one set 
of candidate landforms in Paper 2, but because of geological 
ambiguities it does not reach the 'conclusive' standard. 

If the vapor injection mechanism is a fumarole or gas- 
charged fountain rather than a lake surface, then the ability 
to set a paleopressure constraint goes away. The potentially 
much higher injection rate can overcome the dilution by the 
thicker atmosphere. 

7.2. Applications to the pre-modern sedimentary 
record 

One motivation for our localized-precipitation work is to 
understand Late Hesperian and Amazonian fluvial features 
which are widespread, although rare, on the Martian surface 
(Figure 15; Paper 2). They postdate the global decline in 
valley formation, aqueous mineralisation and erosion rates 
near the Noachian-Hesperian boundary. Therefore, local- 
ized processes are a priori a reasonable explanation for these 
fluvial features. Unusual microclimates (such as nearby 
steep slopes) may explain some of the features [Fassett et 
al, 2010], but those on the Valles Marineris plateau and 
at Mojave Crater suggest localized precipitation (Figure 15; 
Mangold et al. [2008]; Paper 2). 

Impact-induced precipitation has been proposed to ex- 
plain fluvial fans at Mojave Crater, but not yet modeled 
[ Williams & Malin, 2008]. Mojave is a fresh crater whose in- 
ner terraces are dissected by channelized fans that drain hill- 
slopes. Several other Amazonian craters show Mojave-like 
fans, but they are more degraded (e.g., PSP_007447_1995 
in Tartarus and PSP_008669_1905 in Isidis). Mojave has 
been described as a 'Rosetta Stone' for decoding impact- 
associated fluvial features on Mars, because it is relatively 
young and records with high fidelity processes that are de- 
graded in the ancient record [Segura et al, 2008; Senft & 
Stewart, 2008; Toon et al, 2010]. 

7.3. Can localized precipitation explain Noachian 
erosion? 

Early Mars' geomorphic record strongly suggests that 
long-lived and regional/global precipitation contributed to 
Noachian erosion, but it is also consistent with a role for 
localized precipitation. The well-dated Late Noachian / 
Early Hesperian [Fassett & Head, 2008b; Hynek et al., 2010] 
peak in valley formation on Mars includes some integrated 
drainage networks extending over 1000s of km, that appear 
to have formed over a long period of time relative to crater 
formation (e.g., Barnhart et al. [2009]). Because localized 
precipitation on a cold planet is almost always short-lived 
precipitation (§1), these findings are not consistent with 
the model developed in this paper. In addition to this 
timescale problem, the rapid discharges and short formation 
timescales required to form some valley networks and deltas 
are at best marginally compatible with snowmelt [Kraal et 
al, 2008; Kleinhans et al, 2010], and suggest ejecta dewater- 
ing or rain-phase precipitation. (Atmospheric liquid water 
aerosol is currently not included in our model.) 

On the other hand, the majority of drainage systems on 
Early Mars are poorly integrated and localized, including 



most infilled crater floors and eroded crater rims. Inspired 
by observations at Mojave Crater [Williams & Malin, 2008], 
and supported by the results of the idealized model pre- 
sented here, we suggest that localized precipitation is rel- 
evant to understanding these systems. Mars' deuterium- 
hydrogen ratio indicates that Mars has lost water over time 
[Jakosky, 1991]. Therefore, crater lakes resulting from im- 
pacts into icy targets, wet-based ice sheets and lava-ice inter- 
actions would all have been more common on ancient Mars, 
and could serve as localized vapor sources. An Early Mars 
GCM with water sources at the valley-associated crater lakes 
shows precipitation only near the lakes and at planetary 
cold traps, consistent with our mesoscale results [Soto et 
al, 2010]. 

Any hypothesis for Early Mars erosion should explain the 
transition from significant crater infilling and crater rim ero- 
sion prior to the Noachian-Hesperian boundary, to minimal 
erosion in the Mid-Late Hesperian and Amazonian (as doc- 
umented by, for example, Golombek et al. [2006]; Forsberg- 
Taylor et al. [2004]; Boyce & Garbeil [2007]; Howard et al. 
[2005]; Irwin et al. [2005]; Moore & Howard [2005]). In a 
localized precipitation scenario, this would correspond to 
a change in planetary target properties: an ice-rich crust 
before the erosional transition, and an ice-poor crust after- 
wards. 

For impacts into a planetary crust containing buried ice 
layers, models suggest that fast-moving ice-rich debris flows 
cause crater infilling in the minutes after impact [Senft & 
Stewart, 2008]. Therefore, we focus on relatively late-stage, 
relatively prolonged (hours to millenia) erosion of the inner 
slope of the crater rim. 

Energy balance suggests that it is physically possible for 
localized precipitation to contribute to increased Noachian 
crater rim erosion. The energy source for the system is the 
heat of shocked target material. If the vapor source close 
to the crater center (fumarole, geyser or lake) can mine 
heat from deep within the shocked material, a sustained va- 
por source for years is possible, but the interval over which 
snow melts on contact with fresh ejecta on the rim is limited 
by conductive cooling of the surface boundary layer. How- 
ever, if snowmelt-fuelled erosion removes the chilled bound- 
ary layer, hot material will again be exposed and melting 
can resume. Therefore, positive feedback between erosion 
and melting can sustain snowmelt-driven erosional activity 
on the crater rim. Snow can mine heat from deeper within 
the ejecta blanket provided that 



(T b - 273.15)cj > — (L melt + (273.15 - T sn )c sn ) (14) 

where Wt/Rk is the water:rock ratio required for erosion, 
Lmeit is the latent heat of melting, T sn is snow temperature, 
and Can is snow heat capacity. 

Ejecta blanket thickness is variable but lunar scalings sug- 
gest craters > 160 km diameter have ejecta > 0.6 km thick 
at their rim [McGetchin et al, 1973]. Therefore, if ejecta 
temperature exceeds 900K, snow is continuously supplied at 
1.25 mm/hr as found in our simulations, and the majority 
of eroded material is removed by debris flows with a W/R 
of 1.5, condition (14) is satisfied, and > 0.6 km material can 
be eroded from the inner slope of the crater rim by localized 
precipitation within a decade. 

This is directly relevant to three of the final four Mars 
Science Laboratory candidate landing sites. Gale (155 km 
diameter), Holden (155 km diameter), and Eberswalde (65 
km diameter; immediately NE of Holden), because all three 
have deeply dissected rims showing evidence for aqueous 
transport of material from the rim to the crater floor [An- 
derson & Bell, 2010; Moore & Howard, 2005; Rice & Bell, 
2010]. 
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8. Conclusions 

We conclude from this study that:- 

(1) A low-pressure lake effect creates deep convection, 
rapid updrafts, and intense precipitation above cold liquid 
water surfaces on Mars. On Earth, because of much higher 
atmospheric pressure, this requires tropical temperatures. 

(2) We use MRAMS to simulate lake storms on Mars. The 
modeled storms have updraft velocities and plume heights 
which exceed the intensity of the strongest recorded thun- 
derstorms on Earth. 

(3) The fraction of vapor that is trapped near the vapor 
source as localized precipitation increases with lake size. 

(4) The localized greenhouse effect of the released water 
vapor is too weak to cause melting of the snow. 

(5) Melting of equatorial, rapidly-emplaced localized 
snow with the albedo of dust is possible for a subset of or- 
bital conditions, even with the present-day weak greenhouse 
effect. 

(6) Assuming that transient lakes on Mars are uncorre- 
cted with orbital forcing, melting of rapidly-emplaced lo- 
calized precipitation is more likely than melting of precip- 
itation that has been emplaced in equilibrium with orbital 
conditions. 

(7) Taken together, localized storms, rapidly-emplaced 
localized precipitation, and favorable orbital conditions pro- 
vide an alternative working hypothesis for (at least part of) 
the erosion and channel formation observed on pre-modern 
Mars. 

Appendix A: Methods 

The fully compressible dynamical core of MRAMS is de- 
rived from the terrestrial RAMS code. A cloud micro- 
physical scheme derived from the Community Aerosol and 
Radiation Model for Atmospheres (CARMA) was recently 
added to MRAMS by T.I. Michaels. This cloud microphysics 
scheme has been used to successfully reproduce spacecraft 
observations of low-latitude clouds downwind of the Tharsis 
Montes and Olympus Mons [Michaels et al, 2006]. A recent 
description of MRAMS capabilities is Michaels & Rafkin 
[2008b]. 

We made minimal modifications to the MRAMS v. 2. 5-01 
code to allow for surface liquid water. Liquid water micro- 
physics is not included. Water vapor thermodynamics are 
included in the energy equation, but water vapor is not in- 
cluded in the mass and momentum equations - that is, we 
ignore pressure and virtual temperature effects. We do not 
permit dynamic dust lifting at the mesoscale. 

Model vertical layer thicknesses varied from 2.3 km at al- 
titude to 30 m near the ground. We used four grids with 
the outermost being hemispheric and a horizontal resolution 
of ~5.9 km on the inmost grid. Output was sampled every 
1/24 sol (« 3699 s), or 'Mars-hour.' We assume that this 
frequency, limited by available disk space, is enough to cap- 
ture model behavior - for example, we refer to the warmest 
of 24 samples during a sol as the 'day's maximum tempera- 
ture.' The timestep varied between runs but was never more 
than 3.75s on the inmost grid. 

The season is southern summer (Ls ~ 270°) for the runs, 
and boundary conditions are from the NASA Ames MGCM, 
version 2.1. We use present-day orbital conditions. 

Our spinup procedure included 24 Mars-hours with no 
vapor release from the lake, 3 Mars-hours with vapor re- 
lease but cloud microphysics off, and the remainder of the 
run with aerosol microphysics on. We observed that obvious 
spin-up transients had died away by the end of day 3. Snow 
albedo feedback was not considered, but would tend to in- 
crease the intensity of convergence and convective intensity 
by increasing the temperature gradient between land and 
lake. 
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Table 1. List of runs with parameters. 



Run 


Full name 


Lake size (pixels) 


Description 




Deep moist convection? 


dry 


lake_0.005_Mar_7_2010 


- 


Control case without lake 




- 


ref 


lake_0.000_Jan_10_2010 


11x11 


Reference simulation (13 days) 






line 


lake_0.001_Jan_10_2010 


5x149 


Line source in valley 






s 


lake_0.018_Mar_30_2010 


lxl 


Size sensitivity 




X 


m 


lake_0.019_Mar_30_2010 


3x3 


Size sensitivity 




X 


1 


lake_0.030JVlay_10_2010 


7x7 


Size sensitivity 




V 


xxl 


lake_0.020_Mar_30_2010 


29x29 


Size sensitivity 






rough 


lake_0.006_Mar_12_2010 


11x11 


Garrett [1992] lake-surface roug 


;hness 


V 


lo_sun 


lake_0.023_Apr_16_2010 


11x11 


0.75 x present Solar luminosity 




s/ 


hires 


lake_0.004_Feb_10_2010 


Tlxll' 


Horizontal resolution check 




v 7 


winter 


raK.e_u.uuu_ivray l_zuiu 


1 1 y1 1 


oedboiidji aenoiiiviiy \\-is — o\j j 




/ 

V 


midlat 


lake_0.000_Jun_25_2010 


11x11 


Latitude sensitivity (latitude ~ 


45°) 


V 


hipress 


lake_0.033_Jul_22_2010 


11x11 


60 mbar atmospheric pressure 




X 


dz_30m 


afternoon_dz_test_30m_Jan_ll_2010 


11x11 


Vertical resolution check 






dz_10m 


afternoon_dz_test_10m_Jan_ll_2010 


11x11 


Vertical resolution check 






dz_3m 


afternoon_dz_test_3m_Jan_ll_2010 


11x11 


Vertical resolution check 






dz_lm 


afternoon_dz_test_lm_Jan_ll_2010 


11x11 


Vertical resolution check 







Table 2. Vapor fate (end of sol 6). Units are Mt (10 6 metric 
tons). Italicized dry run is subtracted from all other runs. See 
also Figure 6. 



Run 


Evap. rate (mm/hr) 


Water in atm. 


Total atm .(%) 


Snow in lake 


Snow beyond lake 


Total snow (%) 


dry 





420 


100% 








0% 


ref (sol 6) 


2.5 


334 


26% 


234 


741 


74% 


ref (sol 12) 


2.5 


316 


11% 


543 


1942 


89% 


line 


3.7 


1420 


12% 


2701 


7355 


88% 


s 


16.2 


26 


38% 


0.8 


42.7 


62% 


m 


3.6 


68 


49% 


12 


59.1 


51% 


1 


2.5 


228 


44% 


84 


204 


56% 


xxl 


2.9 


899 


9% 


5030 


4503 


91% 


rough 


2.6 


381 


27% 


250 


793 


73% 


lo_sun 


2.8 


176 


12% 


484 


820 


88% 


hipress 


3.4 


175 


22% 


484 


820 


78% 
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Figure 1. To show dependence of lake-driven con- 
vection on lake temperature and atmospheric pressure. 
Moist convection is enhanced on low-pressure Mars rela- 
tive to 1-bar Earth because more buoyancy is produced 
for a given sea-surface temperature or lake-surface tem- 
perature. One measure of the strength of convection is 
the maximum vertical velocity of an updraft w ma x ( = 
\/2CAPE, where CAPE is Convective Available Poten- 
tial Energy), the so-called thermodynamic speed limit 
[Markowksi & Richardson, 2010]. Suppose the Level of 
Free Convection (LFC) to be at the surface and the Equi- 
librium Level (EL) at 1.5km for both planets. Suppose 
complete isobaric condensation and precipitation of a 
parcel at 50% humidity at the surface. In an otherwise 
dry atmosphere, and neglecting some second-order ther- 
modynamic corrections, w ma x is then as shown. The red 
curve is for today's Mars, the blue curve is for today's 
Earth, and the magenta curve is is the pressure used in 
our Mars simulations. The black arrow shows that con- 
vection above a 5°C lake on Mars may be as vigorous as 
above a 72° lake on Earth. The orange curve corresponds 
to the threshhold pressure above which we suspect local- 
ized precipitation on Mars does not occur. Values used: 
c p ,Mars 770 J/kg; c PiEart h 1003 J/kg; R Ea rth = 287 J 
K _1 mol _1 ; Ruars = 189 J K~ x mol -1 ; gMars = 3.7 m 
s~ 2 ; gEarth = 9.8 m s~ 2 ; L sub = L evap ~ 2.5 x 10 6 J/kg; 
saturation vapor pressure curve from [Hardy, 1998]. 
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Ice fraction and total water column 




100 150 200 



Figure 2. Time averaged response to the lake. The 
area shown is a subset of the inmost grid, and distances 
are in km from lake center. The shaded area in the first 
panel shows the extent of the lake, a) Temperature per- 
turbation due to lake; b) radiative forcing at the surface; 
c) overall surface-layer wind field with lake added; b) 
perturbation to surface-layer wind field due to lake; e) 
total water column abundance (thin lines at spacing of 
0.25 precipitable cm) and fractional ice abundance (thick 
lines) . 
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Figure 3. Spatially-averaged time dependencies are 
highly repeatable between sols, (a) Lake perturbation 
to surface temperature downwind of the lake, as a func- 
tion of surface temperature in K. (b) Mass of water in 
atmosphere in 400km side square box centered on the 
cloud, as a function of near-surface air temperature in K. 
Green circles correspond to vapor mass, and black crosses 
correspond to ice mass. 
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65N 75N 85N 

Figure 4. N-S cross section through lake storm. Blue 
tint corresponds to increasing water ice fraction (interval 
0.002, maximum value 0.011). Labelled contours corre- 
spond to bulk vertical velocity in m s _1 . y-axis is vertical 
distance in m. x-axis is horizontal distance in simulation 
units: 10 simulation units = 59 km. Lake extends from 
69N to 79N on this scale. 
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Figure 5. Precipitation from our ref simulation, (a) 
Mean precipitation (mm/hr water equivalent) for refer- 
ence simulation. Peak value is 1.5 mm/hr. (b) Maximum 
precipitation (mm/hr water equivalent) for reference sim- 
ulation. Peak value is 6 mm/hr. (c) Diurnal cycle in 
spatially-averaged precipitation (mm/hr) as a function 
of average surface temperature (K): snow falls at a high 
rate during the night but is reduced following sunrise (di- 
agonal branch), levelling out at values ~ | of the night- 
time peak during the late afternoon (horizontal branch). 
Sunset permits a rapid return to high rates of snowfall 
(vertical branch). 
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Figure 6. Fate of vapor released from the lake (precipi- 
tation efficiency). Colored dots corresponds to a models' 
inventory of water substance after 6 sols, after substract- 
ing the inventory of a lake-free run. Color corresponds 
to fraction of atmospheric water in the ice phase: red 
is more ice-rich, blue is more vapor-rich. Only snow on 
ground can contribute to localized geomorphology. Water 
substance in atmosphere can contribute to regional and 
global climate change: greenhouse warming is increas- 
ingly likely as the mass of atmospheric water increases. 
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Figure 7. Outcome of vertical resolution sensitivity test. 
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Figure 8. Intensity of moist convection for runs listed in 
Table 1. Red diamonds correspond to peak time-averaged 
plume velocity, blue circles correspond to plume height. 
Lines connect the run of 5 simulations in which only 
size was varied. Plume height is defined as the altitude 
at which the plumes' spatial peak in time-averaged ice 
mixing ratio falls below 10 -3 . The two blue circles for 
hipress correspond to ice mixing ratios of 10 -3 (lower 
circle) and 10 -4 (upper circle). The two red diamonds 
for hires correspond to peak velocities at 2km resolution 
(upper diamond) and peak velocities smoothed to 5.9 km 
resolution (lower diamond). 
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Lake area (km 2 ) 



Figure 9. Dependence of plume convective intensity on 
lake size. The black line is a best-fit power law (slope = 
1.11) to the three central points. It intersects the lakeless- 
run Planetary Boundary Layer (PBL) turbulence at a 
lake size of O(10 2 )km 2 ; below this lake size, we would 
not expect to recognize a plume. Explanations of the 
similarity of the power-law exponent to 1, and the devi- 
ation of the largest lake from this slope, are provided in 
the text. 
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Figure 10. Exceedance probabilities for annual max- 
imum sunlight at Mars' equator. Vertical dashed line 
is 640 W/m 2 , the minimum for melting (see text). The 
small wiggles in the solid curves are interpolation arti- 
facts. Snow instantaneously emplaced at Mars' equator 
with the current orbital conditions and solar luminosity 
would melt. For progressively earlier times (reduced solar 
luminosity), the probability of melting decreases. 
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Figure 11. Exceedance probabilities for annual peak 
temperature of snowpack at Mars' equator, from our ID 
column model. A dust-like albedo (0.28) is assumed. Ver- 
tical dashed line corresponds to the melting point. The 
small wiggles in the solid curves are interpolation arti- 
facts. For progressively earlier times (reduced solar lu- 
minosity), the probability of melting decreases. 
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Figure 12. Self-normalized probability density func- 
tions for the subset of sampled orbital conditions that 
produce melt. Albedo = 0.28, age = 1.0 Gya (other albe- 
dos and ages have very similar pdf shapes, although max- 
imum temperatures differ). The cyan lines corresponds 
to 273K (melting), and the other colors are green (all 
temperatures), blue (268K), and red (278K). 
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Figure 13. Contoured probability (percent) of melting 
of rapidly-emplaced equatorial precipitation as a function 
of albedo and solar luminosity. In the case of impact- 
induced precipitation, snow falling on hot ejecta will melt 
regardless of orbital conditions. 
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Figure 14. Seasonal supraglacial lakes on Greenland 
ice sheet. Note absence of clouds and precipitation at 1 
bar, and high drainage density of meandering, snowmelt- 
fed channels feeding lake. Largest lake is 4.5 km across. 
(Image used with kind permission of I. Joughin, U. Wash- 
ington Polar Science Center). 
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Figure 15. Primary motivation for this study: Two 
geological settings where fluvial activity postdates the 
Late Noachian - Early Hesperian acme of channel for- 
mation, and there is a candidate vapor source nearby. 
Scale bars are 500 m. Green box: Sub-kilometer fans re- 
sembling alluvial fans at Mojave Crater, which may be a 
recent impact into icy ground (7.9N 326. 6E, orthorecti- 
fied HiRISE image PSP_001481_1875). Channels dissect 
both sides of 100-200m tall ridge and channel heads are 
found < 20m from ridgeline.; Red box: Inverted streams 
on plateau near Juventae Chasma (4.3S 296. 7E, orthorec- 
tified HiRISE image PSP_003223_1755). Negligible post- 
Noachian erosion over most of the rest of the planet 
indicates that such events were localized. Sources for 
stratigraphic context: Hartmann [2005] (absolute ages); 
Murchie et al. [2009] (clay and sulfate stratigraphy) Carr 
& Head [2010] (age of outflow channels); Fassett & Head 
[2008b] (age of valley networks); Masse et al [2010] & 
Mangold et al. [2010] (existence of young sulfates); Le 
Deit et al. [2010] (age of Juventae plateau channels). 



